Universal Coefficient of Performance at Maximum Figure of Merit for Tight-Coupling 

Refrigerators 
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We find and prove that for a refrigerator tightly and symmetrically coupled with two heat baths 
at temperatures T c and Th (> T c ), the coefficient of performance at maximum figure of merit 
asymptotically approaches to y/ec when the relative temperature difference between two heat baths 
Eq 1 = (Th — Tc)/T c is small enough. This universal coefficient of performance at maximum figure 
of merit can be regarded as the counterpart of universal efficiency at maximum power output for 
tight-coupling heat engines operating between two heat baths at small temperature difference in the 
presence of left-right symmetry. 

PACS numbers: 05.70.Ln 



Introduction-One of key topics in finite-time thermo- 
dynamics is the efficiency at maximum power (EMP) for 
heat engines. Since the pioneer work made by Curzon 
and Ahlborn this topic has been fully investigated by 
many researchers [2h20||. Recent researches are mainly 
focused on two key issues: One is the universal EMP for 
tight-coupling heat engines operating between two baths 
at small temperature difference [t]|4l3|: the other is the 
global bounds of EMP for heat engines operating between 
two baths with arbitrary temperature difference [13j-|20| . 
In particular, Van den Broeck [fj found that the universal 
EMP for tight-coupling heat engines was equal to ?yc/2 
up to the first order term of relative temperature differ- 
ence between two baths, where r\c is the Carnot efficiency 
which can also be understood as the relative temperature 
difference for heat engines. The universality up to the 
second order term of rjc was first observed in Ref. Q, 
and then proposed as a conjecture in Ref. and finally 
proved by Esposito et al. for tight-coupling heat engines 
in the presence of left- right symmetry (lpj . In addition, 
Esposito et al. also found the EMP of heat engines to 
be bounded between rin/2 and ?yc/(2 — r)c) under low- 
dissipation conditions [13]. It is interesting that these 
two bounds are also shared by the linear-irreversible heat 
engines [lj|. The accessibility of the bounds for differ- 
ent kinds of heat engines are also fully investigated 0- 
[2fjj . These researches highly enhance our understanding 
of the applicable scope of the Curzon- Ahlborn efficiency 
t]ca = 1 — V 1 ~~ Vc f°r endoreversible heat engines. 

On the other hand, the optimal performance of refrig- 
erators has also attracted much attention 21 -33|. Since 
the power input is not an appropriate optimal target 
function for Carnot-like refrigerators, it is relatively dif- 
ficult to define an optimal criterion and obtain its cor- 
responding coefficient of performance (COP) for refrig- 
erators as researchers did for heat engines. Velasco et 
al. [1^] adopted the per-unit-time COP as a target func- 
tion and proved Eca = V £ c + 1 — 1 to be the upper 
bound of COP for endoreversible refrigerators operat- 
ing at the maximum per-unit-time COP, where Ec is 



the Carnot COP for reversible refrigerators. Yan and 
Chen [30] suggested to take sQc/t C ycic as the target func- 
tion, where t cyc \ c is the time for completing the whole 
Carnot-like cycle, e and Q c are respectively the COP 
of refrigerators and the heat absorbed by the working 
substance from the cold bath. They also optimized this 
target function and found the corresponding COP to be 
Eca = V e c + 1 — 1 f° r endoreversible refrigerators [3(| . 
Recently, de Tomas et al. introduced a unified optimiza- 
tion ^-criterion for heat devices which was defined as 
the product of the conversion efficiency and the heat ab- 
sorbed per unit time by the working substance [3l| . This 
X-criterion returns to the power output for heat engines 
while degenerates into the criterion proposed by Yan and 
Chen for refrigerators. Additionally, the COP at maxi- 
mum x was also found to be eca = VI + £ c — 1 for 
symmetric low-dissipation refrigerators [3l|. Based on 
the work by de Tomas et al., one of the present authors 
and his coworkers derived that the COP at maximum x 
was bounded between and (y/9 + 8ec — 3)/2 for low- 
dissipation refrigerators [32|. The observed COP's for 
real refrigerators are also located in the region between 
these bounds, which is in good agreement with their the- 
oretical estimation. These bounds were also confirmed 
by Izumida et al. with a minimally nonlinear irreversible 
model for refrigerators 



The above researches on refrigerators support that x~ 
criterion is an appropriate figure of merit for refrigera- 
tors. The results obtained from maximizing this figure 
of merit for refrigerators have counterparts in those de- 
rived from maximizing the power output for heat engines. 
To begin with, the COP at maximum figure of merit 
(eca = \/l + £c~ 1) for the Yan-Chen endoreversible re- 



vT+eo-1) 
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for the Curzon- Ahlborn endoreversible heat engines [30| . 
Next, the bounds of COP [0 and (y/9 + 8e c - 3)/2] at 
maximum figure of merit for low-dissipation refrigera- 
tors correspond to the bounds of efficiency [?yc/2 and 
ryc/(2 — tic) ] at maximum power for low-dissipation heat 
engines |32J. In particular, the COP at maximum fig- 
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ure of merit (sca = VI + £ c — 1) for symmetric low- 
dissipation refrigerators also corresponds to the EMP 
(tjca = 1 — VI — Vc) fo r symmetric low-dissipation heat 
engines j3f[. However, it is still lack of the counterpart 
of universal EMP for tight-coupling heat engines in the 
COP at maximum figure of merit for tight-coupling re- 
frigerators when they operate between two baths at small 
temperature difference. Our main goal in this work is 
to complement this theoretical imperfection. We opti- 
mize the figure of merit for the Feynman ratchet as a 
typical tight-coupling refrigerator operating between two 
baths and then find that its COP at maximum figure of 
merit approaches to ^fec (Note that e^ 1 can also be un- 



derstood as the relative temperature difference) when it 
operates between two baths at small temperature differ- 
ence. It is obvious to see that this behavior is also shared 
by £ca = V 1 + £ c — 1) which implies that y/e~c might 
be the universal COP at maximum figure of merit for 
tight-coupling refrigerators working between two baths at 
small temperature difference. This proposition is proved 
to be true under certain symmetric conditions within 
the framework of linear irreversible thermodynamics by 
adopting the new convention on the thermodynamic flux 
related to the heat transfer between two baths. 

Feynman ratchet as a refrigerator.-The Feynman 
ratchet can be simplified as a Brownian particle 
walking in a periodic lattice labeled by © n , (n — 
• • • , —2,-1, 0, 1, 2, • • • ) with the fixed step size 8. The 
ratchet potential is schematically depicted in Fig. [TJ 
where the energy scale and the position of potential 
barrier are respectively denoted by e and 56. The pa- 
rameter S, so called load distribution factor, takes value 
in between and 1. The Brownian particle is in con- 
tact with a cold bath at temperature T c in the left 
side of each potential barrier while it in contact with 
a hot bath at temperature Th (> T c ) in the right side 
of each barrier. The particle is pulled from the left to 
the right side by a moment z due to the external force. 
In steady state, the forward and backward jumping rates 
can be respectively expressed as w+ = koe~( e ~ zSe ^ Tc and 
uj- = koe~( e ~ zSe+ze ^ Th according to the Arrhenius law 
34j . In the expressions of jumping rates, the Boltzmann 
factor is set to 1 while fco represents the bare rate constant 
with dimension of time . For simplicity, we introduce 
two abbreviated notations q = e — z59 and w = z8. 

For the relative large load z, the forward jumping rate 
can be larger than the backward one. In this case, the 
net current 



J 



UJ + 



-q/T c 



-(q+w)/T h 



(1) 



is positive. In each forward step, the particle absorbs 
heat q = e — z58 from the cold bath. Combining the in- 
put work w = zO done by the external load, the absorbed 
heat will be released into the hot bath when the particle 
jumps over the barrier. Thus the total heat q + w will be 




FIG. 1. Schematic digram of Feynman ratchet as a refriegra- 
tor. 



released into the hot bath in each forward step. The en- 
ergy conversion in each backward step is exactly opposite 
of that in forward step mentioned above. Thus the net 
power input can be expressed as P — w J while the heat 
absorbed from the cold bath or released into the hot bath 
per unit time can be expressed Q c — qj or Qh = (q+w) J, 
respectively. Obviously, when J > the heat flows from 
the cold bath to the hot one, and the power input (i.e., 
the mechanical flux) is proportional to the thermal fluxes 
(Q c and Qt). It is in this sense that the Feynman ratchet 
is regarded as a tight-coupling refrigerator. 

The COP of this tight-coupling refrigerator can be ex- 
pressed as 



Q c /P = q/w. 



(2) 



Simultaneously, the figure of merit can be expressed as 
k q 2 



X = eQc 



-q/T c 



-(q+w)/T h 



(3) 



Maximizing x with respect to the internal barrier hight 
e and the external load z, we can obtain 

(2 - q/T c )e-*/ T ° = (2- q/T h )e-W*, (4) 
e-« /T c = (I + w/T h )e- {q+w),Th . (5) 

Combining Eq. ^ and the above two equations, we de- 
rive that the COP at maximum figure of merit (e*) sat- 
isfies the following transcendental equation: 

E£jz±(i-±) = ^;^c (6) 

where e c = T c /(T h - T c ) is the COP of Carnot refrig- 
erators. It is interesting and surprising that depends 
merely on ec (or equivalently speaking, the relative tem- 
perature difference between two baths) rather than the 
load distribution factor 5 although the expressions of 
jumping rates contain 5. 

It is almost impossible to achieve the analytic solution 
to Eq. ([6]). We will investigate the asymptotic behaviors 
at large temperature difference (ec — > 0) and small tem- 
perature difference (ec —> oo), respectively. For the for- 
mer case, considering < e* < ec, we transform Eq. (J5]) 
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into 2ec/e* + £*/ec — 3 = ln(2ec/£*), from which we 
obtain 

e, ~ 0.524e c (7) 

when — >• 0. This behavior is different from that of 
£ca = V £ c + 1 — 1 which gives eqa — 0.5ec when 
ec —> 0. Thus there is no universal behavior of COP at 
maximum figure of merit for refrigerators at large tem- 
perature difference. 

On the other hand, it is not hard to prove £* — > oo and 
£*/sc ~^ when the temperature difference between two 
baths is very small (i.e., ec — > oo) [35]. If multiplying 
1 + e^ 1 on both sides of Eq. © and then expanding it 
into a series of e~ x and £q, we can derive 

2/e c - £*/4 + 2/e*£c - 2/el = 0, (8) 

when we neglect the contribution of higher order terms. 
The terms £*/e c an d 2/e*ec are of higher order relative 
to 1/sc since e» — > oo and £*/ec — > when ec —> oo, 
thus Eq. ((5J is further simplified into 2/ec — 2/e* — 
when ec — > oo. Its solution is 

e* ~ y/ec~, (9) 

which gives the asymptotic behavior of COP at maximum 
figure of merit for the Feynman ratchet as a refrigerator 
operating between two heat baths at small temperature 
difference. 

We suggest to use an interpolation formula 

e* = \J ec + a 2 — a (10) 

with a = 1/(2 x 0.524) = 0.954 as the approximate so- 
lution to Eq. ©. It is easy to see that this formula 
degenerates into Eq. ([7]) and Eq. ([9]) when ec — > and 
oo, respectively. We compare this interpolation formula 
with the numerical solution to Eq. © in Fig. [2] Sur- 
prisingly, this interpolation formula (solid line) does ex- 
tremely approach to the numerical solutions (squares) to 
Eq. ([5]) obtained from the high precision computation 
when ec takes values in a relatively large range. This 
formula fits the numerical data better than formulas ([7]) 
and © depicted respectively as dot line and dash dot 
line in Fig. There exists a little difference between 
this interpolation formula and ecA = V £ c + 1 — 1 ■ Both 
of them are respectively depicted as the solid line and 
dash line in the inset of Fig. [21 which reveals that the 
interpolation formula is much closer to the numerical so- 
lutions than ecA = V £ c + 1 — 1. The small relative 
error (<0.8%) between interpolation formula (flU)) and 
the genuinely exact solution to Eq. is quantitatively 
estimated in Ref. [35| . 

Universality. -It is obvious that ecA = V £ c + 1 — 1 — 
yfec when ec — > oo. That is, the Yan-Chen endorc- 
versible refrigerators or the symmetric low-dissipation re- 
frigerators also share the same asymptotic behavior © 
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FIG. 2. (Color online) COP at maximum figure of merit for 
the Feynman ratchet as a refrigerator. The numerical solu- 
tions to Eq. ([6]) and interpolation formula (|10[) are depicted 
as the squares and solid line, respectively. The diagrams of 
functions (0 and (|9} are depicted as the dot line and dash 
dot line, respectively. Inset graph shows the diagram of inter- 
polation formula (solid line) and that of sca = Vsc + 1 — 1 
(dash line) as well as the numerical solutions (squares) in the 
range of 5 < ec < 24. 



as the Feynman ratchet at small temperature difference. 
Since all of them can be regarded as tight-coupling refrig- 
erators, we may conjecture that a universal COP at max- 
imum figure of merit e* ~ ^fec exists for tight-coupling 
refrigerators working between two baths at small tem- 
perature difference. 

As it was done by Van den Broeck [f| for heat engines, 
here we consider a generic setup for a tight-coupling re- 
frigerator shown in Fig. [3l An external force F is ap- 
plied on the system and inputs a power P = Fx into 
the system, where x is the thermodynamically conjugate 
variable of F. The dot represents the derivative with re- 
spect to time. The corresponding thermodynamic force 
is X\ — F/T, where T is the temperature of the system 
which can be well defined due to the assumption of lo- 
cal equilibrium. The thermodynamic flux conjugated to 
X\ is Ji = x. Then the power input can be expressed 
as P = J\X\T in terms of the thermodynamic flux and 
force. Assume that the system is in contact with a cold 
bath at temperature T c = T — s c AT and a hot bath at 
temperature T h = T+s h AT with AT < T. The positive 
parameters s c and sh should satisfy s c + sh — 1 due to 
the constraint AT = Th — T c . Their specific values de- 
pend respectively on the coupling strengths between the 
model system and the cold or hot baths. 

In steady state, the entropy production rate can be ex- 
pressed as a — Qh/Th — Q c /T c which can be further ex- 
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pressed as a = (Q c + P)/T h -Q c /T c = P/T h +Q c (l/T h - 
1/T C ) or o = Q h /T h - (Q fc - P)/T c = P/T f + Q h (l/T h - 
1/T C ) due to the conservation of energy (Q c + P = Q h ). 
Both equations can be respectively expressed as a — 
(P/T)[l - s h AT/T + 0(AT/T) 2 } +Q c (l/T h - 1/T C ) and 
a = {P/T)[l + s c AT/T + 0{AT/T) 2 }+Q h (l/T h -l/T c ) 
with the consideration of Th = T + s^AT and T c 
7 - s c AT, where 0(AT/T) 2 represents the term in the 
same order of (AT/T) 2 . Since s c + Sh = 1, we can fur- 
ther derive ct = s c cr + S/l cr = (P/T)[l + C(AT/T) 2 ] + 
(s c Qc + ShQh)(^-/Th — l/?c) which enlightens us to take 
X 2 = 1/T h - 1/T C and 

J% = s c Q c + shQh (if) 

as the thermodynamic force and flux related to the heat 
transfer between two heat baths, respectively. The trun- 
cation error with the consideration of new convention 
([TTj) is in the order of (AT/T) 2 when the entropy pro- 
duction rate is expressed as a = J\X\ + J2X2. The 
other conventions such as J2 = Q c or J2 = Qh result 
in the lower accuracy with the truncation error in the 
order of AT/T when the entropy production rate is ex- 
pressed as a = J\X\ + J 2^-2- Therefore, we will adopt 
the high-precision convention to deal with the tight- 
coupling refrigerator. With the consideration of Eq. (fTTj) 
and Q c = Qh — P, the heat absorbed from the cold bath 
can be expressed as Q c = J<i — ShP ■ 

? T h 

Qh 

^\^^ ^ 

Qc 

1 Tc 

FIG. 3. Generic setup for a tight-coupling refrigerator. 

According to linear irreversible thermodynamics, we 
write the linear relationship between the thermodynamic 
fluxes and forces: 

J\ = L\\X\ + J2 = L/2\X\ + -^22-^2, (12) 

where the Onsager coefficients satisfy L\\ > 0, L22 > 0, 
L11L22 — -^12-^21 > and L\2 = L^i- Furthermore, the 
tight-coupling condition L\ 2 — L\ x = L11L22 leads to 
J2/J1 = Lu/Lu. With the consideration of this equa- 
tion and P = TJ\Xi, the COP can be expressed as 
e = Qc/P = L 2 i/TLnXi — Sh, from which we have 
X\ = L\2/TL\\(e + Sh). Substituting this equation into 



the definition of figure of merit, we obtain 

X = eQ c - L 22 e 2 [l + TX 2 (e + s h )]/T(e + s h ) 2 . (13) 

Maximizing \ with respect to e (equivalently to X\), 
we find that the COP at maximum figure of merit sat- 
isfies el + 3s h e* + 2s 2 h + 2sh/TX 2 — which gives 
e* = ^s 2 /4 - 2s h /TX 2 - 3s h /2. With the considera- 
tion of X 2 = 1/T h - 1/T C and e c = T c /(T h - T c ), we 
achieve 

£* = y/e%/4 + 2s h e c (e C + l)/(e c + s c ) - 3s h /2. (14) 

If the model system is symmetrically coupled with both 
heat baths such that Sh = s c = 1/2, we derive e* ~ y/sc 
when sc — > 00 from Eq. (|14p . Therefore, we have proved 
the conjecture on the universal COP at maximum figure 
of merit for tight-coupling refrigerators symmetrically in- 
teracting with two heat baths at small temperature dif- 
ference. In addition, in the extremely asymmetric cases 
of Sh — (s c = 1) and Sh — 1 (s c = 0), equation (fT4")) 
leads to £o = and E\ = (-y/9 + 8ec — 3)/2, respectively, 
which surprisingly equal to the global lower and upper 
bounds [e_ = and e+ = (V9 + 8e c - 3)/2] of COP 
at maximum figure of merit for low-dissipation refriger- 
ators which are also found to be reached in the case of 
extremely asymmetric dissipations pjil ]. 

Conclusion and discussion.- In the above discussion, 
we investigate the COP at maximum figure of merit for 
the Feynman ratchet as a refrigerator and find that the 
corresponding COP can be approximately expressed as 
interpolation formula (|10j) . In the limit of small tem- 
perature difference between two baths, this formula and 
£ca = v 7 ! + £ c — 1 share the same asymptotic behav- 
ior (jnj). This universal asymptotic behavior is proved to 
be available for refrigerators tightly and symmetrically 
coupled with two baths at small temperature difference. 

The key point in our proof of universality is to adopt 
new convention (|lll) as the thermodynamic flux related 
to the heat transfer between two baths. If we adopt this 
convention and recalculate the model for a tight-coupling 
heat engine in Ref. jg], we still obtain the same result of 
universal EMP (rjc/2) as Van den Broeck did. Interest- 
ingly, if further considering the symmetric-coupling be- 
tween the engine and its two baths, we can return to 
the universal EMP [r]c/2 + r] c /8 proved by Esposito et 
al. up to the second order term 35]. In this sense, 

asymptotic behavior ^ for refrigerators tightly and sym- 
metrically coupled with two baths at small temperature 
difference can be regarded as the counterpart of universal 
EMP for tight-coupling heat engines operating between 
two baths at small temperature difference in the presence 
of left-right symmetry. In addition, we find the same sur- 
prise that the global lower and upper bounds of EMP for 
low-dissipation heat engines obtained in Ref. [13] can be 
achieved when we adopt convention (fTTj) to deal with the 
tight-coupling heat engine in the extremely asymmetric 
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cases of s c — (si = 1) and s c = 1 (s/j = 0) [35| , It is a 
fantastic fact that a local theory valid for the small tem- 
perature difference can give the correctly global bounds 
available for the arbitrary temperature difference. In our 
opinion, it is a challenge to investigate whether this fact 
takes place merely by coincidence or due to some under- 
lying reasons. 
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